Simulation of Quantum Adiabatic Search in the Presence of Noise 
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Results are presented of a large-scale simulation of the quantum adiabatic search (QuAdS) al- 
gorithm in the presence of noise. The algorithm is applied to the iVP-Complete problem N-Bit 
Exact Cover 3 (EC3). The noise is assumed to Zeeman-couple to the qubits and its effects on the 
algorithm's performance is studied for various levels of noise power, and for 4 different types of 
noise polarization. We examine the scaling relation between the number of bits N (EC3 problem- 
size) and the algorithm's noise- averaged median run-time (T(N)). Clear evidence is found of the 
algorithm's sensitivity to noise. Two fits to the simulation results were done: (1) power-law scaling 
(T(N)) = aN b ; and (2) exponential scaling (T(N)) = a [exp(foiV) - 1]. Both types of scaling rela- 
tions provided excellent fits, although the scaling parameters a and b varied with noise power, and 
with the type of noise polarization. The sensitivity of the scaling exponent b to noise polarization 
allows a relative assessment of which noise polarizations are most problematic for quantum adiabatic 
search. We demonstrate how the noise leads to decoherence in QuAdS, and estimate the amount 
of decoherence present in our simulations. An upper bound is also derived for the noise- averaged 
QuAdS success probability in the limit of weak noise that is appropriate for our simulations. 

PACS numbers: 03.67.Lx,03.65.Yz,05.40.Ca 
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I. INTRODUCTION 



In computational complexity theory [jj , computational 
problems are classified according to the resources needed 
to obtain a solution. Often such a problem is reformu- 
lated as a decision problem whose solution is a "yes" or 
"no" answer. An algorithm that solves all instances of 
a decision problem is said to be a polynomial-time algo- 
rithm if the time needed to find a solution to an arbitrary 
instance is 0(n k ), where n is the size of the instance 
and k is a fixed positive integer. The computational 
complexity class P is composed of decision problems for 
which polynomial-time algorithms exist. Another im- 
portant computational complexity class is the class NP 
which is composed of decision problems for which, for 
each such problem, a polynomial-time algorithm exists 
to verify the "yes" output of a candidate solution. The 
question of whether the classes P and NP are equal is 
one of the biggest open problems in theoretical computer 
science. It is widely conjectured that P ^ NP. A deci- 
sion problem p is said to be polynomially transformable 
to a decision problem q if: (1) there exists a function / 
that maps every instance x of p into an instance f(x) 
of q for which the answer to x is "yes" if and only if 
the answer to f(x) is "yes"; and (2) a polynomial-time 
algorithm exists to compute f(x) for every x in p. A 
problem is said to be iVP-Complete if it belongs to NP 
and every problem in NP is polynomially transformable 
to it. Thus, if a polynomial-time algorithm is found for 
an iVP-Complete problem, it follows that P = NP. In 
light of Shor's polynomial-time quantum factoring algo- 
rithm 0, Bennett et. al. asked whether a polynomial- 
time quantum algorithm might exist to solve an NP- 
Complete problem. In other words, they asked, "Can a 
quantum computer solve an ./VP-Complete problem in 
polynomial-time?" Should the answer to this question 



prove to be "yes" , it is widely recognized that this would 
have profound consequences for theoretical computer sci- 
ence Q|. 

In 2001, Far hi et. al. Q examined the question of 
whether a quantum computer might be able to solve an 
./VP-Complete problem in polynomial-time. They ap- 
plied the quantum adiabatic search (QuAdS) algorithm 
introduced in Ref. to find solutions of randomly gen- 
erated hard instances of the ./VP-Complete problem N- 
Bit Exact Cover 3 which they believe to be classically 
intractable for sufficiently large inputs. The algorithm 
was simulated numerically on an existing (classical) com- 
puter. Because a quantum computer cannot be simulated 
efficiently on a classical computer, the simulations were 
restricted to 7 < N < 20. They found that the median 
run-time T(N) for QuAdS to succeed on this class of in- 
stances could be fit with a quadratic scaling relation in 
the number of bits N (problem-size): T(N) ~ N 2 . It was 
pointed out that should classical algorithms truly require 
exponential time on this class of instances, and should 
the quadratic scaling behavior of QuAdS persist to large 
N, then QUAdS could outperform classical algorithms 
on randomly generated hard instances, though not nec- 
essarily on worst case instances. Their work has drawn a 
great deal of interest and suggests two possible directions 
for further research: (1) examining the large N scaling 
behavior of QuAdS; and (2) examining the robustness of 
QuAdS to noise. Recent progress on (1) appears in the 
work of Roland and Cerf @; Mitchell et. al. @; and La- 
torre and Oriis 0, ITol ] . In this paper we will focus on the 
question of robustness of QuAdS to noise. 

The analysis of Ref. p| assumed that the quantum 
computer was perfectly isolated from noise. Childs et. al. 
[Tl| were the first to consider the robustness of QuAdS to 
noise. They presented numerical simulation results which 
suggested QuAdS might have an inherent robustness to 
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two types of noise- induced errors: (1) transitions out of 
the instantaneous ground state; and (2) unitary control 
errors. Their simulation did not, however, attempt to 
match the scale of the simulation presented in Ref. 
Recently, Roland and Cerf carried out an analytical study 
of QuAdS in the presence of noise using random matrix 
theory They showed that the error probability of 

QuAdS would not increase with problem-size N if the 
noise varies either very slowly or very rapidly with re- 
spect to the natural time-scale of the quantum computa- 
tion h/E, where E is the energy scale for the eigenvalues 
of the instantaneous Hamiltonian H(t). Their analysis 
assumed weak noise which justified a perturbative anal- 
ysis of their model. 

As stated above, in this paper we will also examine 
the robustness of QuAdS to noise. Our work comple- 
ments the analysis of Ref. ^1 m that we consider noise 
that varies on a time-scale comparable to Ti[E. We ex- 
tend the simulation protocol of Farhi et. al. @ to include 
noise, and determine the noise-averaged median runtime 
(T(N)} for QuAdS to succeed on randomly generated in- 
stances of N-Bit Exact Cover 3 which have a unique solu- 
tion. We find clear evidence of the sensitivity of QuAdS 
to noise. Two fits to the simulation results were carried 
out: (1) power-law scaling (T(N)) = aN b ; and (2) ex- 
ponential scaling (T(N)) = a [exp(WV) — 1]. Both types 
of scaling relations provided excellent fits, although the 
scaling parameters a and b varied with noise power, and 
with the noise polarization (defined in Section Ull Aj l. Al- 
though we did encounter examples of noise realizations 
that reduce the runtime for QuAdS to succeed, at the 
largest noise power levels that we examined, the pre- 
dominant effect of noise was to slowdown QuAdS. The 
sensitivity of QuAdS to noise polarization allows, for the 
first time, a determination of which noise polarizations 
are most problematic for QuAdS. 

The structure of this paper is as follows. In section ITTl 
we briefly describe the ./VP-Complete problem N-Bit Ex- 
act Cover 3, and review how QuAdS can be used to find 
a solution to an instance of this problem. In section ITTT1 
we describe our noise model and explain our extension of 
the simulation protocol of Farhi et. al. to include noise. 
In section IIVI we present the results of our simulation, 
and discuss their significance in section [3 Section ]V\ 
also demonstrates how noise-induced dephasing leads to 
decoherence in QuAdS, estimates the amount of deco- 
herence present in our simulations, and derives an upper 
bound for the noise-averaged QuAdS success probability 
in the weak noise limit that is appropriate for our simu- 
lations. Finally, in section IVT1 we provide a summary of 
our work. 



II. EXACT COVER 3 AND QUANTUM 
ADIABATIC SEARCH 

In this section we: (i) present the A^P-Complete prob- 
lem Exact Cover 3 (section III : (ii) describe the quan- 



tum adiabatic search algorithm and how it is used to 
solve instances of Exact Cover 3 (section IIIBfl ; and (iii) 
describe the simulation protocol of Ref. [j| fsection fll Cfl . 

A. Exact Cover 3 

Consider a collection of N binary variables (bits) 
zi, ■ ■ ■ , zn each of which can take on the value or 1. 
The state of this AT-bit system is specified by assigning 
values to each bit in the string z = z% ■ ■ ■ Zjj. A total of 
2 N states z are possible. 

A specific instance of Exact Cover 3 (EC3) is composed 
of L clauses, each of which imposes a constraint on the 
values of three of the N bits. The number of clauses L 
will generally vary from one EC3-instance to another. If 
the i-th clause involves the bits (a(i), b(i), c(i)), then a 
particular bit-string z = z% - ■ • Zjy satisfies the clause if 
z a(%) + z b{i) + z c{%) = 1- Otherwise, one says that z does 
not satisfied the clause. A bit-string z is said to solve an 
instance of EC3 if it satisfies all the clauses that make 
up the instance. If one is given an instance of EC3 and 
a candidate solution string z* , one can check whether 
z* solves the instance in polynomial-time since checking 
each of the L clauses simply requires adding three inte- 
gers, and this can be done in polynomial-time. Thus EC3 
belongs to NP. It can be shown that all problems in NP 
are polynomially transformable to EC3 [jjj so that EC3 
is also AP-Complete. 

B. QuAdS and EC3 

By combining the dynamics of the quantum adiabatic 
theorem with a clever choice of Hamiltonian H(t), quan- 
tum adiabatic search (QuAdS) causes the state of a 
quantum computer (QC) to home-in on a solution of 
an instance of an ./VP-Complete problem We as- 
sume the AP-Complete problem is EC3. The quan- 
tum adiabatic theorem [lj| assures us that if the initial 
state of a quantum system is the groundstate \E g (0)) 
of the initial Hamiltonian H(0), and H(t) varies suffi- 
ciently slowly, then the quantum state at time T can 
be brought arbitrarily close to the groundstate of H(T): 
\tp(T)) = \E g (T)) + \5ip). The probability not to be in the 
groundstate, e(T) = \{5il)\5ijj)\ 2 , vanishes as T — > oo. For 
a given instance of EC3, QuAdS evolves the state \tp(t)} 
over a time interval < t < T, and the final Hamiltonian 
H(T) — H p is constructed so that its groundstate en- 
codes a solution of the given instance. The initial Hamil- 
tonian H(0) = En is chosen so that its groundstate is 
non-degenerate and can be easily constructed, and H(t) 
smoothly morphs Hi — ► H p : 

ff(*)=(l-|) H i+Cf) H v ■ (!) 

We now briefly describe how Hi and H p are constructed 
for a given instance of EC3. 
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Suppose that we are given an instance of A-bit EC3 
composed of L clauses Cj = (a(i), b(i), c{i) ) with 
i = 1, • • ■ , L, and that the bit-string s — si ■ ■ ■ sn is a 
solution of this instance. Our QC is assumed to contain 
N qubits and the computational basis states \z\ - ■ • Zjf) 
are chosen to be eigenstates of o\ ® • • • ® . Although 
the QuAdS algorithm can be formulated more abstractly, 
it proves convenient to adopt the language of NMR. 
Hi, then, describes the Zeeman coupling of the qubits 
to an external magnetic field which points along the x- 
direction. The strength of the magnetic field at the site 
of the i-th qubit is equal to the number of clauses in 
the EC3 instance that contain bit i. The non-degenerate 
groundstate of Hi thus has all qubit spins aligned along 
x so that \E g (0)) is simply the uniform superposition of 
all 2^ computational basis states \zx ■ ■ ■ zm)- The final 
Hamiltonian H p is constructed to be diagonal in the com- 
putational basis \z\ ■ ■ ■ zn)'- 



the measurement result z — s. The algorithm's failure 
probability e(T) will be much less than 1 if |14j : 



H p \z%--- z N ) = h(zi, • ■ • , z N ) \zi ■ ■ ■ z N ) 



(2) 



The eigenvalue h(zi, • ■ • , zjv) is the sum of energy func- 
tions h Ci (z a (i), z b (i), z c{i) ), 

L 

h{zi,--- ,zn) = /] hCi(z a (i), Zb(j), z c (j)) , (3) 



where C, = (a(i), b(i), c(i)) is the i-th clause in the given 
EC3 instance and, 



h-C t {Za(i),Z b (i),Z c (i)) 



if zx • ■ • zn satisfes Cj 

1 otherwise . 



Thus h(zx, ■ • • , zn) indicates how many clauses are vio- 
lated by the string z = z\ ■ • ■ Zpf. All strings s = Si ■ ■ -sjv 
that solve the given instance correspond to computa- 
tional basis states si ■ • ■ sjv) that have zero energy, and 
which together span the groundstate eigenspace of H p . 
EC3 instances which have a unique solution are referred 
to as unique satisfying assignment (USA) instances, and 
are believed to be the most difficult for QuAdS. In 
essence, it's harder to find a needle in a haystack that 
contains only one needle, than it is to find one in a 
haystack that contains many needles. For further dis- 
cussion of instances with multiple satisfying assignments, 
see Ref. [H|. 

Thus, if our QC is initially prepared in the groundstate 
of Hi, and T is chosen to be sufficiently large that H(t) 
evolves adiabatically, then the final state \ip(T)) will lie 
in the groundstate eigenspace of H p with probability 1 — 
e(T), where e(T) = \{6ip\5ip}\ 2 < 1. For a USA instance 
with solution string s = si • • • sjsr, measuring |-0(T)) in 
the computational basis will give the solution string s 
with probability 1 — e(T). Let the string z\ ■ ■ ■ zn be the 
actual measurement result. It can be quickly tested to 
determine whether it solves the USA instance. If "yes" , 
the algorithm has succeeded, and z = s. If "no", then 
the quantum adiabatic search procedure is repeated until 



where, 



e 
A 



max 

t e [0,T] 



T > — 
A 2 



( El (t)\T^-\E g (t)) 



t wn \Exit) - E g (t)\ ; 



(4) 

(5) 
(6) 



and \E g (t)) (\Ei(t))) is the groundstate (first excited 
state) of H(t). 



C. Noiseless Simulation Protocol 

Farhi and co-workers 0, 0] studied the performance 
of QuAdS by numerically integrating the Schrodinger 
equation using H(t) (eqn. Q) to drive the dynamics. 
Because the Hilbert space for Af-qubits has dimension 
2 N , numerical simulation quickly becomes impractical 
since the number of wavefunction components that must 
be tracked by the numerical integration grows exponen- 
tially with N. This practical difficulty caused the simu- 
lations done in Refs. 0, 0| to be restricted to the range 
7 < N < 20. Only randomly generated USA instances 
were simulated as these authors believed that this set 
would provide hard cases for QuAdS, though not neces- 
sarily the worst cases. 

An A-bit USA instance of EC3 is generated by the fol- 
lowing procedure. The first clause is generated by picking 
3 distinct integers in the range [1, N] at random (uni- 
form deviate). A count of the number of bit-strings of 
length N that satisfy the clause is then done. Distinct 
clauses continue to be generated in this manner, and with 
each new clause generated, a new count of the number 
of bit-strings that satisfy all clauses generated to that 
point is done. Eventually, the point is reached where 
enough clauses have been generated so that one of two 
situations occurs. (1) Only one bit-string s of length N 
remains that satisfies all the clauses generated. The col- 
lection of clauses then corresponds to an EC3-instance 
which has a unique solution s, and the procedure has 
thus produced a USA instance of EC3. (2) Addition of 
the most recently generated clause causes the set of all 
clauses produced to have no satisfying assignment (i. e. 
an EC3-instance has been produced with no solution) . In 
this case, the instance is discarded and the above proce- 
dure is repeated until a USA instance is generated. The 
USA instance produced is then used to construct Hi and 
Hp, and from them, H(t). The USA instances generated 
by this procedure were found to contain approximately 
as many clauses L as number of bits N: L ~ N . 

In the adiabatic limit, T — > oo, the QuAdS success 
probability P s — > 1. Numerical simulations must nec- 
essarily work with finite T so that a protocol is needed 
to determine how to pick T. Farhi et. al. carried out 
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FIG. 1: Simulation results generated by our QuAdS numeri- 
cal code for the noiseless median runtime T(N) (dimensionless 
units) versus the number of bits N. The solid line is the best 
power-law fit to the data and the dash-dot line is the best 
exponential fit. Respective values for a and b are given in 
the text. The error bars give 95% confidence limits for each 
median. Each datapoint is the outcome of averaging over 75 
USA instances. 



a hunting procedure which searches for a T that causes 
the simulation to produce a success probability P s in the 
range [0.12, 0.13]. This corresponds to a success proba- 
bility in the vicinity of P s = 0.125 = l/2 :i . As noted in 
Ref. [j| , this value of P s is somewhat arbitrary, although 
for the values of N used in the simulations, it is much 
larger than 1/2 N which is the probability of the solution 
string s in the initial state. For each value of N, 75 USA 
instances are generated and QuAdS is used to find the 
solution to each instance. This generates runtimes Ti(N) 
{i = 1, • • • , 75) from which the median runtime T(N) is 
determined. Farhi et. al. found that their simulation 
results for T(N) could be fit with a quadratic scaling re- 
lation for the range of A^- values considered: T(N) ~ A^ 2 . 

As a check on the soundness of our own simulation 
code, we repeated the noiseless calculation of Farhi et. al. 
for 7 < N < 14. Figure [I] shows our results for T(N) 
versus N, together with a power-law fit to the simula- 
tion data: T(N) = aN b . The best- fit parameters are 
a = 0.11966 and b = 2.0034. The x 2 value for the fit is 

0. 092, and the probability P(v 2 > 0.092) = 0.9999844. 
As noted in Section llVl the closer this probability is to 

1, the more consistent the data-set is with the fitting 
function. We see that a power-law fit provides excellent 
agreement with our simulation results, and our exponent 
b = 2.0034 is consistent with the quadratic fit of Ref. 5]. 
Figure n a l so includes an exponential fit to the noiseless 
data. The best-fit parameters are: a = 3.89566; b = 
0.140235; x)u = 0- 214 ; and p (x 2 > 0.214) = 0.9998121. 
We see that exponential scaling also provides an excellent 
fit to the data. We will comment further on this in Sec- 
tions IIVI and [V] The addition of noise to the simulation 
of QuAdS causes simulation to become impractical when 
N > 12. In a couple of cases we were able to simulate 



N = 13 and 14 (see Section Hvjl . Our noiseless simulation 
will be used as a baseline for our simulation of QuAdS in 
the presence of noise. Since the simulations with noise do 
not go beyond N = 14, we did not extend the noiseless 
simulations beyond N = 14. 



III. NOISY QUANTUM ADIABATIC SEARCH 

In this section we discuss QuAdS in the presence 
of noise. Section IIII Al presents our noise model and 
describes how noise realizations are constructed. Sec- 
tion ^^B] discusses how the noiseless simulation protocol 
is modified to allow a study of the noise-averaged perfor- 
mance of QuAdS. 



A. Noise Model 

We consider noise Ni(t) that couples to the qubits via 
the Zccman interaction: 

N 

i=i 

The qubits are labeled by the index i and the subscript 
i on the noise field Ni(i) and coupling constant 7* indi- 
cates that they can vary from one qubit site to another. 
We assume the noise is a stationary random process. To 
simplify the presentation of our noise model we initially 
consider a one-component noise field N(t). The presen- 
tation is easily generalized to three-component spatially 
varying noise Nj(f). 

The rate at which the noise field N(t) can do work 
(i. e. noise power) is |l6j. 

P = N 2 (t) , 

and the energy that can be delivered in a time interval 
dt is, 

dE = N 2 {t)dt . 

We consider power-type noise for which the time- 
averaged noise power 

P= lim i f T/2 N 2 (t)dt (8) 

is finite. The total noise energy 

/oo 
dtN 2 {t) (9) 
-00 

diverges for this class of noise. The divergence is due to 
the occurrence of an infinite number of noise fluctuations 
in the time interval — 00 < t < 00. The energy of an 
individual fluctuation is, however, finite. 
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The mean noise power P can be related to the noise 
correlation function, 

1 f T ' 2 

N(t)N(t~ s) ee lim - / dyN(y)N(y-s) . (10) 



-T/2 



Comparing eqns. (JHJ and (|1C)|> we see that, 



P = N 2 (t) 



(11) 



The Weiner-Khintchine theorem [ljj shows that the 
noise correlation function and the power spectral density 
SV(/) form a Fourier transform pair: 



N(t)N(t-s) = dfS N (f)e- 2 ^ s . (12) 

J —OG 

Thus, it follows from eqns. (jllfl and (|12[) that 

/* / 4f<M/) , (13) 



which identifies Sn(/) as the mean noise power available 
in the frequency interval (/, / + df). 

Having reviewed these basic facts about stationary ran- 
dom processes, we are now in a position to describe our 
noise model and to write out the basic relations that are 
necessary to construct an individual realization of the 
noise. 

The noise N(t) is produced by a sequence of randomly 
occurring noise fluctuations F(t). The fluctuations: (1) 
occur independently of each other at average rate n per 
unit time; and (2) have a peak value x which is Gaussian 
distributed with mean x = 0, variance x 2 — a 2 , and 
temporal width 2t, where r is the thermal relaxation 
time. The bandwidth of F(t) is Alu ~ l/2r. Thus, 



N(t) = ^f(i- U 



where i labels the noise fluctuations, and £j specifies the 
center of the zth fluctuation. The mean number of fluctu- 
ations Aff occurring in the time interval [0, T] is Aff = nT. 
It is well-known that for noise with these properties, the 
actual number of fluctuations n that occur in a time T is 
governed by the Poisson distribution 0] : 



The energy present in a single fluctuation is: 

F 2 {t)dt . 



(14) 



Let F(t) — xh(t), where h(t) is any convenient function 
of finite support with normalization 



dth 2 {t) = 2t 



As mentioned above, x is Gaussian distributed with mean 
x = and x 2 = a 2 . From eqn. Ijl4|l . £ = 2x 2 r, and the 
mean energy per fluctuation e is, 

e = 2x^T = 2a 2 T . (15) 

From Campbell's theorem [l^, the power spectral den- 
sity for N(t) is 



S N (f) = n\g(f)f 



(16) 



where <?(/) is the Fourier transform of the fluctuation 
profile F(t). Thus, using eqns. (IS} , (HH) , and Paresval's 
theorem gives, 



(17) 



P = n dtF 2 {t) 



Finally, using eqns. (|14|l and 115(1 gives, 

P = 2na 2 r . (18) 

Thus we see that our noise model is characterized by any 
3 of the parameters P, n, a 2 , and r. 

The numerical simulation constructs a realization of 
the noise as follows. We sample a positive integer Nf ac- 
cording to the Poisson distribution with mean Aff — nT, 
where T is the duration of the search. Nf is the num- 
ber of fluctuations present in the noise realization. The 
noise model assumes these fluctuations occur indepen- 
dently with probability dpf = (1/T)dt. We sample Nf 
numbers t% (i = 1, • • ■ , Nf) from the interval (0, T). The 
U correspond to the temporal centers of the Nf fluctua- 
tions. For simplicity, we've assumed that the fluctuation 
profile h(t) is a square pulse. We next carry out Nf sam- 
ples Xi (i = 1, • • • , Nf) of a Gaussian distribution with 
mean xi — and variance x 2 = a 2 . Here Xi is the peak 
value of the i-th fluctuation. These sample results pro- 
duce the noise realization N(t): 



N(t) 



N f 

x * 



sgn(i - tu) - sgn(< - t ir ) 



(19) 



where tu = ti — T, and U r = U + t. We shall need to pro- 
duce noise realizations with arbitrary mean noise power 
P. We do this by the following normalization procedure. 
First we calculate the mean noise power V of the noise 
realization N(t) just produced: 



V 



1 

f 



dtN 2 (t) 



(20) 



Then, if the desired value for the noise power is P, 
we rescale N(t) in eqn. (|19f) so that N(t) — > n(t) = 



P/V N(t). The result is a noise realization n(t) with 
mean noise power P. The simulation takes as inputs 
the mean noise power P, the variance x 2 = a 2 , and the 
thermal relaxation time r. The fluctuation rate n then 
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follows from eqn. 1|18|) : n = P/(2a 2 r). This procedure 
is used to produce 4 types of noise that are character- 
ized by the polarization, or direction along which the 
noise Ni(i) fluctuates: (i) x; (ii) y; (iii) z; and (iv) all 3 
directions simultaneously. We shall refer to these noise 
polarization types as x-type; y-type; z-type; and 3-type 
noise, respectively. The noise Ni(i) is then introduced 
into eqn. J7J and the full Hamiltonian TL(t) is the sum of 
H(t) (eqn. JTJ}) and H mt {t) (eqn. ©): 

H(t)=H(t)+H mt (t) . (21) 

H(t) drives the Schrodinger dynamics of QuAdS in the 
presence of the noise realization Ni(i). We now go on to 
explain how the noiseless simulation protocol of Ref. |5j 
is extended to include the effects of noise. 



B. Noisy Simulation Protocol 

We would like to compare how the 4 different noise 
types introduced above impact the performance of 
QuAdS. To that end, we determine the scaling relation 
for the noise-avcragcd median runtime (T(A)} versus N 
for each type of noise. Just as in the noiseless proto- 
col, for each value of N, we generate 75 USA instances of 
EC3. QuAdS, in the presence of each of the 4 noise types, 
is then applied to the same 75 USA instances. This al- 
lows an apples-to-apples comparison of how the different 
noise types affect QuAdS performance. For each noise 
type, and each USA instance, we generate 10 noise re- 
alizations and use the Hamiltonian TL{t) to find the 10 
corresponding QuAdS runtimes T^{N). Here j labels 
the noise realizations (j = 1, ■•■ ,10), and i labels the 
USA instances (i = 1, • • • , 75). Thus, for each value of 
N, and each noise type, 750 runtimes are found. The 
ensemble En of 750 runtimes T^(N) allows us to esti- 
mate the statistical effects of each noise type on QuAdS 
performance over the full set of 75 USA instances. The 
noise-averaged median runtime (T(N)) is then identified 
with the median runtime calculated from the noise en- 
semble £m- We then fit the simulation results using both 
power- law scaling (T(N)) — aN b and exponential scal- 
ing (T(N)) = a [exp(WV) _ i] j an d calculated the x 2 for 
each type of fit. We determine 2 such fits for each of the 
4 noise types. For each fit-type, we compared the respec- 
tive scaling exponents to get a quantitative measure of 
which of the 4 noise types most adversely affects QuAdS. 
As will be discussed in Section IIVI scaling curves will 
be determined for each of the 4 noise types for 3 dif- 
ferent noise power levels. As pointed out at the end 
of Section 111 CI most simulations will be restricted to 
7 < N < 12, though for x-typc noise, the upper limit 
could be extended to N = 13 and 14. It is important 
to recognize that each point on a given scaling curve is 
based on 750 runtimes, and so each scaling curve with 
7 < N < 12 (7 < N < 14) is distilled from 4500 (6000) 
runtimes. Each runtime T({N) is itself the result of a 



hunting procedure that requires, on average, 5 integra- 
tions of the Schrodinger equation. Thus a single scaling 
curve corresponds to approximately 22,500 (30,000) in- 
tegrations of the Schrodinger equation. As mentioned 
above, we generate 4 such curves for each noise power 
level, and we present data for 3 power levels. Putting 
all this together, we see that the simulation results we 
present in this paper are the outcome of approximately 
300,000 integrations of the Schrodinger equation. This 
work was done, initially on our own 16-node Beowulf 
cluster, and later on the TeraGrid cluster which was ac- 
cessed through the National Center for Supercomputing 
Applications (NCSA) in Urbana, Illinois. 



IV. SIMULATION RESULTS 

In this Section we present our simulation results for 
noisy QuAdS. The noise model parameters (see Sec- 
tion IIII|I were chosen to be a = 0.2; r = 1; and 
P = 0.001, 0.003, 0.005. This produced noise that was 
strong enough to affect QuAdS performance, though not 
so strong as to make large-scale simulation impractical. 
Both Hi and H p have energy-level spacing AE ~ 1 . For 
our choice of t, AE ~ 2/r, so that our noise bandwidth 
matches the natural resonance energies of the quantum 
computer. We are thus probing QuAdS in a noise regime 
that complements that studied in Ref. Further dis- 
cussion of our choice of noise parameters will be given in 
Section [3 As this was our first large-scale simulation of 
noisy QuAdS, we decided to simplify the noise model by 
restricting the noise field and coupling constants to be 
the same for each qubit: Ni(i) = N(i) and 7^ = 1. Wc 
plan to lift this restriction in our next set of simulations. 
We also plan to parallelize the hunting procedure which 
should substantially speed up the code, allowing us to ex- 
plore larger values of mean noise power P. To get a sense 
of how long the following simulation took in real time, 
note that producing one scaling curve at P = 0.005 took 
approximately 3 weeks on our 16-node Beowulf cluster. 
The same simulation took approximately 72-96 hours on 
the TeraGrid cluster at NCSA. 

We now present our simulation results. For each mean 
noise power P we present 4 figures and 2 tables. Each of 
the figures contains our numerical results for (T(N)) for a 
particular noise polarization type (x, y, z, 3), along with 
two fits to the data: (i) a power-law fit (T(N)) ~ aN b ; 
and (ii) an exponential fit (T(N)) — a [exp(6A) — 1]. 
Information about the best fits to the data are collected 
in the two tables, one table for each of the fitting func- 
tions. Each table contains: (i) the best fit parameters; 
(ii) the chi-squared for the fit xjiv an< ^ ( m ) the prob- 
ability P(x 2 > Xfit) that, assuming the fitting function 
correctly describes the scaling of (T(N)) versus N, a sam- 
pling of (T(N)) would yield a \ 2 > X%t- Tne closer this 
probability is to 1, the more consistent the data-set is 
with the fitting function. As such, it provides a quan- 
titative measure of the quality of fit to the data. Note 



that, due to the finite-size of the data-set, more than one 
fitting function can be consistent with the data. To help 
the reader locate specific results, we provide the following 
roadmap through the figures and tables. 



1. Figures 13 El 0| and [5] correspond to x, y, z, and 

3-type noise with P = 0.001, respectively. TablesH] 
and [n] give the specifics of the power-law and ex- 
ponential fits through these data-sets, respectively. 
Table entries are ordered according to increasing 
scaling-exponent b. 



2. Figures El IH1 and [3] correspond to x, y, z, and 3- 
type noise with P = 0.003, respectively. Tables ITTT1 
and IIVI give the specifics of the power-law and ex- 
ponential fits through these data-sets, respectively. 
Table entries are ordered according to increasing 
scaling-exponent b. 



3. Figurcs HUlimir^l audlT^lcorresDond to x. y, z, and 

3-type noise with P = 0.005, respectively. Tables IVl 
and IVII give the specifics of the power-law and ex- 
ponential fits through these data-sets, respectively. 
Table entries are ordered according to increasing 
scaling-exponent b. 



Examination of these figures and tables shows that the 
power-law and exponential fits both provide excellent fits 
to the numerical results. One sees that the probability 
P{x 2 > X%t)' an d hence the quality of fit, shows a slight 
decrease with increasing mean noise power P, and that 
the rate of decrease in the fit quality is largest for y-type 
noise. One also notes that the quality of the exponential 
fit decreases more slowly than does the quality of the 
power-law fit. Further discussion of these results is given 
m Section |V| 



TABLE I: Summary of best-fit parameters for power-law scal- 
ing curves for P = 0.001. For comparison, best-fit parameters 
for noiseless QuAdS are also included. 



p 


noise type 


a 


b 


Xfit 


P(X 2 > Xfit) 


0.000 




1.1966 x 10" 1 


2.0034 


0.092 


0.9999844 


0.001 


X 


1.1447 x 10" 1 


2.0267 


0.083 


0.9999886 


0.001 


3 


1.0635 x lO^ 1 


2.0594 


0.064 


0.9994915 


0.001 


z 


1.0395 x 10 _1 


2.0695 


0.063 


0.9995085 


0.001 


y 


9.6766 x 10~ 2 


2.1075 


0.078 


0.9992617 




N 

FIG. 2: Simulation results for QuAdS with noise polarized 
along x and mean noise power P = 0.001. Plotted are the 
noise-averaged median runtime {T(N)} (dimensionless units) 
versus the number of bits N. The solid line is the best power- 
law fit to the data and the dash-dot line is the best exponential 
fit to the data. The error bars give 95% confidence limits for 
each median. Each datapoint is the outcome of averaging over 
75 USA instances and 10 noise realizations per USA instance. 




N 

FIG. 3: Simulation results for QuAdS with noise polarized 
along y and mean noise power P = 0.001. Plotted are the 
noise-averaged median runtime {T(N)) (dimensionless units) 
versus the number of bits N. The solid line is the best power- 
law fit to the data and the dash-dot line is the best exponential 
fit to the data. The error bars give 95% confidence limits for 
each median. Each datapoint is the outcome of averaging over 
75 USA instances and 10 noise realizations per USA instance. 



V. DISCUSSION 

A. Summary 

We examine here what can be learned from the re- 
sults presented in Section IIVI To begin, we note that 
both power-law scaling (T(N)) — aN b and exponential 
scaling (T(N)) = a [exp(6iV) — 1] provide excellent fits 
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— Power-Law Fit 




N 

FIG. 4: Simulation results for QuAdS with noise polarized 
along z and mean noise power P = 0.001. Plotted are the 
noise-averaged median runtime {T(N)) (dimensionless units) 
versus the number of bits N. The solid line is the best power- 
law fit to the data and the dash-dot line is the best exponential 
fit to the data. The error bars give 95% confidence limits for 
each median. Each dataDoint is the outcome of averaging over 




N 

FIG. 5: Simulation results for QuAdS with noise polariza- 
tion along all 3 directions and mean noise power P = 0.001. 
Plotted are the noise-averaged median runtime (T(N)) (di- 
mensionless units) versus the number of bits N. The solid 
line is the best power-law fit to the data and the dash-dot 
line is the best exponential fit to the data. The error bars 
give 95% confidence limits for each median. Each datapoint 
is the outcome of averaging over 75 USA instances and 10 
noise realizations per USA instance. 



to our simulation results, although the fit parameters are 
noise dependent. Tables HI I VII show that the scaling expo- 
nent b increases with increasing mean noise power P, and 
the rate at which it increases depends on the noise po- 
larization, or direction along which the noise field Ni(i) 
fluctuates. Clearly, as the scaling exponent increases, 
QuAdS performance decreases non-linearly. We simu- 
lated noise with noise polarization: (i) along x; (ii) along 



TABLE II: Summary of best-fit parameters for exponential 
scaling curves for P — 0.001. 



p 


noise type 


a 


b 




P(x 2 > Xm) 


0.001 


X 


3.81403 


0.142561 


0.216 


0.9998065 


0.001 


3 


2.60771 


0.172331 


0.075 


0.9993111 


0.001 


z 


2.46414 


0.176985 


0.059 


0.9995798 


0.001 


y 


2.42677 


0.179578 


0.085 


0.9991179 




N 



FIG. 6: Simulation results for QuAdS with noise polarized 
along x and mean noise power P — 0.003. Plotted are the 
noise-averaged median runtime {T(N)) (dimensionless units) 
versus the number of bits N. The solid line is the best power- 
law fit to the data and the dash-dot line is the best exponential 
fit to the data. The error bars give 95% confidence limits for 
each median. Each datapoint is the outcome of averaging over 
75 USA instances and 10 noise realizations per USA instance. 

y; (iii) along z; and (iv) along all 3 directions simul- 
taneously. At the noise power levels considered, x-type 
noise was found to have the smallest impact on QuAdS 
performance, while y-type noise was found to cause the 
largest slowdown. As noted in Section Hvl although the 
quality of both the power-law and exponential fits is ex- 
cellent in all cases considered, we did see a slight re- 
duction in the quality of fit as noise power increased for 
all noise-types except x-type noise. We also noted (see 
Tables |V] and I VII especially) that the quality of the ex- 
ponential fit decreased at a slightly slower rate than did 
the power-law fit. It would clearly be of interest to ex- 
tend the simulations to larger noise power levels to ex- 
amine: (i) how rapidly the scaling exponent b increases 
with noise power; and (ii) how quickly the quality of fit 
for both power-law and exponential fits deteriorates with 
increasing noise power. Given that a quantum computer 
will eventually crossover to classical behavior with suffi- 
cient noise power, and should classical algorithms truly 
require exponential time on the randomly generated hard 
instances considered in Ref. [f| and here, then one would 
anticipate a power-law fit to eventually become inconsis- 




N 

FIG. 7: Simulation results for QuAdS with noise polarized 
along y and mean noise power P = 0.003. Plotted are the 
noise-averaged median runtime (T(N)) (dimensionless units) 
versus the number of bits N. The solid line is the best power- 
law fit to the data and the dash-dot line is the best exponential 
fit to the data. The error bars give 95% confidence limits for 
each median. Each datapoint is the outcome of averaging over 
75 USA instances and 10 noise realizations per USA instance. 




N 

FIG. 8: Simulation results for QuAdS with noise polarized 
along z and mean noise power P = 0.003. Plotted are the 
noise-averaged median runtime {T(N)) (dimensionless units) 
versus the number of bits N. The solid line is the best power- 
law fit to the data and the dash-dot line is the best exponential 
fit to the data. The error bars give 95% confidence limits for 
each median. Each datapoint is the outcome of averaging over 
75 USA instances and 10 noise realizations per USA instance. 



tent with the simulation data at sufficient noise power 
and the exponential fit to provide the better fit. Beyond 
this threshold value of noise power, QuAdS would then 
be expected to have performance comparable with clas- 
sical algorithms. Observation of such a crossover would 
be extremely important as it would give a direct measure 
of how much noise QuAdS can tolerate. In the following 
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N 

FIG. 9: Simulation results for QuAdS with noise polariza- 
tion along all 3 directions and mean noise power P = 0.003. 
Plotted are the noise-averaged median runtime (T(N)) (di- 
mensionless units) versus the number of bits N. The solid 
line is the best power-law fit to the data and the dash-dot 
line is the best exponential fit to the data. The error bars 
give 95% confidence limits for each median. Each datapoint 
is the outcome of averaging over 75 USA instances and 10 
noise realizations per USA instance. 



TABLE III: Summary of best-fit parameters for power-law 
scaling curves for P = 0.003. For comparison, best-fit param- 
eters for noiseless QuAdS are also included. 



p 


noise type 


a 


b 


V 2 
Xfit 


P(X 2 > Xfit) 


0.000 




1.1966 x lO -1 


2.0034 


0.092 


0.9999844 


0.003 


X 


9.3604 x 10 -2 


2.1304 


0.052 


0.9999971 


0.003 


z 


8.4209 x 10~ 2 


2.1774 


0.050 


0.9996975 


0.003 


3 


7.8825 x 10" 2 


2.2077 


0.083 


0.9991666 


0.003 


y 


2.9821 x 10" 2 


2.6792 


0.570 


0.9663310 



subsection we show that noise-induced dcphasing leads 
to decoherence in QuAdS, estimate how much decoher- 
ence is present in our simulations, and derive an upper 
bound for the noise-averaged QuAdS success probability. 

B. Noise-Induced Decoherence in QuAdS 

Although our focus in this paper is on QuAdS, the 
following analysis can be adapted to the more general 
situation of noisy quantum adiabatic dynamics. We will 
report on this elsewhere. 

The total Hamiltonian for noisy QuAdS is (see 
cq. <EU), 

H(t) = H(t)+H int (t) , 

with H(t) and Hi nt (t) given in eqs. JTJ and Q, respec- 
tively. In the Schrodinger picture, the dynamics is driven 
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TABLE IV: Summary of best-fit parameters for exponential 
scaling curves for P — 0.003. 



p 


noise type 


a 


b 


Xfit 


P(x 2 > xm) 


0.003 


X 


3.24309 


0.157716 


0.169 


0.9999056 


0.003 


z 


2.10817 


0.193125 


0.036 


0.9998428 


0.003 


3 


2.05131 


0.195818 


0.075 


0.9993165 


0.003 


y 


0.847514 


0.285049 


0.283 


0.9908794 




FIG. 10: Simulation results for QuAdS with noise polarized 
along x and mean noise power P = 0.005. Plotted are the 
noise-averaged median runtime {T(N)) (dimensionless units) 
versus the number of bits N. The solid line is the best power- 
law fit to the data and the dash-dot line is the best exponential 
fit to the data. The error bars give 95% confidence limits for 
each median. Each datapoint is the outcome of averaging over 
75 USA instances and 10 noise realizations per USA instance. 



by 



Transforming to an interaction-like picture (the time- 
ordering symbol is suppressed), 



\ip(t)) = exp 



i I dsH mt {s) 
lo 



\x(t)) 



the equation of motion for |x(i)) is found to be 

i^\ x ) = H{t)\ x ) , 

where 

H(t) = e^o dsH int (s) H ^ e -if*dsH int ( 3 ) 
Using the well-known identity 

e^ A B e"« A = B + £ [A, B] + L [A, [A, B]\ + ■ 



(22) 



(23) 



(24) 



(25) 



39 - 
36 - 

33 - 
30 - 
27 - 
A24- 

z, 
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V 

18 - 
15 - 
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3, 



Power-Law Fit 
Exponential Fit 



FIG. 11: Simulation results for QuAdS with noise polarized 
along y and mean noise power P = 0.005. Plotted are the 
noise-averaged median runtime {T(N)) (dimensionless units) 
versus the number of bits N. The solid line is the best power- 
law fit to the data and the dash-dot line is the best exponential 
fit to the data. The error bars give 95% confidence limits for 
each median. Each datapoint is the outcome of averaging over 
75 USA instances and 10 noise realizations per USA instance. 




FIG. 12: Simulation results for QuAdS with noise polarized 
along z and mean noise power P — 0.005. Plotted are the 
noise-averaged median runtime {T(N)) (dimensionless units) 
versus the number of bits N. The solid line is the best power- 
law fit to the data and the dash-dot line is the best exponential 
fit to the data. The error bars give 95% confidence limits for 
each median. Each datapoint is the outcome of averaging over 
75 USA instances and 10 noise realizations per USA instance. 



with B = H(t) and 

£A = if dsH int (s) 
Jo 



— «7 / x(s) ds 
Jo 



(26) 
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FIG. 13: Simulation results for QuAdS with noise polariza- 
tion along all 3 directions and mean noise power P — 0.005. 
Plotted are the noise-averaged median runtime (T(N)) (di- 
mensionless units) versus the number of bits N. The solid 
line is the best power-law fit to the data and the dash-dot 
line is the best exponential fit to the data. The error bars 
give 95% confidence limits for each median. Each datapoint 
is the outcome of averaging over 75 USA instances and 10 
noise realizations per USA instance. 



TABLE V: Summary of best-fit parameters for power-law 
scaling curves for P = 0.005. For comparison, best-fit pa- 
rameters for noiseless QuAdS are also included. 



p 


noise type 


a 


b 


Xfit 


P(X 2 > Xfn) 


0.000 




1.1966 x 10" 1 


2.0034 


0.092 


0.9999844 


0.005 


X 


6.679 x 10 -2 


2.2901 


0.132 


0.9996794 


0.005 


z 


3.5862 x 10" 2 


2.5964 


0.389 


0.9833500 


0.005 


3 


2.4823 x 10~ 2 


2.7752 


0.484 


0.9750332 


0.005 


y 


6.6069 x 10~ 3 


3.429 


1.052 


0.9017628 



identifies 



and 



£ = —ij / x(s) ds 
Jo 



N 



(27) 



(28) 



To make our discussion concrete, I will assume the pres- 
ence of x-noise: Nj(£) = x(t)x; and set ji = 7. Combin- 
ing eqs. lt2^|) through (|2*S|) gives 



H{t) = H(t) -[if J x(s)ds 



N 



,i=l 



In the simplest approximation, only the first term in 
eq. (|29|l is kept. This corresponds to the limit of weak 



TABLE VI: Summary of best-fit parameters for exponential 
scaling curves for P — 0.005. 



p 


noise type 


a 


b 


Xfit 


P(X 2 > Xfit) 


0.005 


X 


2.25018 


0.189744 


0.178 


0.9993283 


0.005 


z 


1.12768 


0.258381 


0.256 


0.9924477 


0.005 


3 


0.859322 


0.287436 


0.253 


0.9926288 


0.005 


y 


0.498037 


0.357153 


0.707 


0.9504644 



noise which is appropriate for our simulations. Clearly 
cq. i|29[) lends itself to a systematic evaluation of cor- 
rections to the weak noise approximation. Evaluating 
cq. lf2*2"|) at the end of the quantum search t = T gives 



IVCH) = exp 



1 I dsH int (s) 
la 



\X(T)) , (30) 



and in the weak noise approximation, eq. I|23|l becomes 
. d 



1-^\X) = H{i)\x) 



(31) 



Thus, in the weak noise limit, the time-development of 
|x(i)) is driven by the Hamiltonian H{t) for noiseless 
QuAdS. 

To evaluate eq. (|3U|) . we introduce the instantaneous 
eigenvalues and eigenstates of H int (t) : 



H int (t)\E a (t)) = E„(t)\E„(t)) 



(32) 



Following the usual decoherence language |2jj, we refer 
to the \Ea-(t)) as the pointer basis states. For x-noise, 
the pointer basis states are time- independent: 



\E a (t)) = \E„ 

= ki 
= M 



• <?n)x 



\<?n)x 



(33) 



where <Tj = ±1 and er = (<7i, . . . , <tjv)- It follows from our 
expression for Hi nt (t) (eqs. JjJ, |3*2*|l . and |3*3*jl ) that 



E a {t) = -"jx{t)o , 



with 



i=l 



(34) 



(35) 



Dividing the time interval [0, T] into N subintervals of 
duration e = T/N, and noting that [H int (t), H int (t')] = 
for x-noise allows us to write the exponential in eq. Ij^Ojl 

as 



exp 



1 I dsH int (s) 
la 



N 



Y[u(u,ti-i) 



(36) 
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where 



?7(*i,*i-i) = exp 



i I dsH int (s) 
hi-i 



(37) 



Using eq. I|36|) and inserting the completeness relation for 
the pointer basis states {l-E^)} into eq. I|30[) gives 



N 



Using eqs. (|3*7)l and (f3*2")l gives 

(<Ti • ■ - CTjv|C/"(ii, = 



(err • • -a N \ exp 



-i / dsE„(s) 



(39) 



Finally, using eq. in eq. gives 
|V(T))= £ |a 1 --- ( x A r}expH^(7 1 )]a<x(r) , (40) 



noise. We now show that these stochastic phases lead to a 
suppression of the off-diagonal matrix elements (<x ^ <x') 
of p(T) when the noise-average is carried out. The noise- 
averaged density matrix p CT CT , (T) is thus 



Pa 



,{T) = D rT ^a -{T)a* rT ,{T) 



where 



= exp [-iT^,] 



(45) 



(46) 



is the decoherence factor. In the adiabatic limit, the 
thermal relaxation time (a.k.a. noise correlation time) r 
satisfies t <SC.T . We divide the integration interval [0, T] 
appearing in T a j (eq. (|44[l) into M — T/t subintervals 
of duration r. This renders r CT CT / into a sum of uncorre- 
cted random variables r CT CT '(j): 



M 

E 



<(j) , 



(47) 



where 



where 



<t>a{T) = [ dsE„(s) = -7(7 / x{s)ds , (41) 



and 



ao-{T) = (o-i • • -(J N \x(T)) 



(42) 



The final density matrix p(T) — \ip(T)){ip(T)\ has matrix 
elements (in the pointer basis) 



Pa 



• (T) = exp [-iT Vt „.] atr (T)aUT), (43) 



where 



r , = mt)-4>At) 



-j (<7 — a') j x(s)ds 
Jo 



(44) 



As noted above, in the weak noise approximation, 
eq. (|31|) indicates that a CT (T) is determined by the noise- 
less QuAdS dynamics. Thus if \ip(0)) = x(0)) 1S initially 
equal to the initial groundstate \E g (0)), then \x(T)) = 
\E g (T)) and so a„{T) = {<j x ■ ■ ■ a N \E g (T)). Thus all 
noise dependence in this approximation appears in the 
stochastic phases {(/>„■ (T)} which clearly depend on er 
(see eq. J410 ). The stochastic character of the noise re- 
quires us to represent our quantum system by an ensem- 
ble in which each element of the ensemble is our quantum 
system in the presence of a particular noise realization. 
As we do not know which element of the ensemble will 
correspond to our quantum system on a given run of 
QuAdS, we must average over the ensemble to determine 
the expected performance of QuAdS in the presence of 



T (T>cr ,(j) = - 1 (a-a') f T x{s)ds . (48) 

•Ai-l)r 

Since the noise is stationary, the set of {T CT CT '(i)} have 
identical probability distributions. If the {r CT)(T '(j)} are 
not only uncorrelated, but also statistically independent, 
it follows from the central limit theorem that r^^./ will 
have a Gaussian probability distribution with mean V a <y i 
and variance r 2 a ,. From eqs. I)47|) and (|48|l . and the fact 

that the noise x(t) has zero mean x{t) = 0, it follows that 
r CT CT ' = 0. From eq. 



T pT 

ds / ds'x(s)x(s') 
li Jo 



Averaging over the noise gives 



7 2 (cr-a'Y 



ds / ds'x(s)x(s') 
ii Jo 



(49) 



(50) 



For our noise model, the noise correlation function 



x(s)x(s') is (see Section ITlI A|) : 



x(s)x(s') = <7 2 h(s — s') 



(51) 



To avoid confusion with a defined in eq. Ij35(l . we have 
written a 2 , for the variance of x-noise (denoted by <7 2 
in Sections IIIII and llV|) ; and recall that h(s — s') is our 
square pulse noise fluctuation profile of unit height and 
width 2t. Using eq. (|51J) in eq. (|50|l and carrying out the 
integrations gives: 



r 2 , 

CT.CT 



4tT {a-a' f^al 



(52) 
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We now have the ingredients needed to evaluate the de- 



coherence factor D CT CT / : 



A, 



exp 



27iT* 



2r 2 , 



exp 



Using eq. fity gives 



exp \-iY a ^i\ 
(53) 

(54) 



A number of remarks are in order. 

(1) For quantum adiabatic dynamics, T > 1. Thus for 
non-zero er 2 , r, and 7, the decoherence factor -D CTiCT ' is 
exponentially small. Thus a quantum state that is a 
superposition of pointer basis states will undergo an 
effective wavefunction collapse and p(T) will become 
effectively diagonal in the pointer basis. Thus we get 
the standard decoherence phenomenology |2£j from our 
noise model. As we have seen, the noisy background 
field Nj(f) generates a stochastic phase <f>o-(T) which is 
different for the different pointer basis states \E„-{t)), 
and which causes suppression of the off-diagonal matrix 
elements of p(T) when we average over the noise. 
As noted above, since one doesn't know what noise 
realization will appear on any particular run of QuAdS, 
one must noise-average to determine the impact of 
noise on QuAdS performance. It is worth noting that, 
although the original decoherence analysis considers 
an environment which is a quantum system, it has been 
conventional wisdom that noise-induced dephasing can 
generate the usual decoherence phenomenology |2lj . 
although we are not aware of any explicit demonstration 
of this expectation prior to the one presented above. 

(2) For QuAdS, the above analysis is relevant whenever 
[Hi n t (t), H(t)] 0, since then the pointer basis states dif- 
fer from the instantaneous eigenstates of H(t). Since the 
instantaneous eigenstates must then be a superposition of 
the pointer basis states, they are vulnerable to the effec- 
tive wavefunction collapse produced by the noise-induced 
decoherence discussed above. Since such an effective col- 
lapse signals the effective loss of quantum coherence, one 
would expect that this decoherence would jeopardize any 
advantage QuAdS might acquire due to the quantum na- 
ture of its dynamics. It was this naive expectation that 
motivated us to do a careful study of the effects of noise 
of QuAdS performance. To test this idea, we decided 
to study noise which coupled to the qubits via the Zee- 
man interaction (see eqn. 0), and whose polarization 
was constant: Nj(t) = x; y; z. The pointer basis states 
are then, respectively, the eigenstates of direct products 
of the qubit Pauli operators: Yl i=1 €3 cr* ; Yl^ =1 ® u % y ; and 



Yli—i ® <j\- We also recognized that such an interaction 
would also allow a study of what would happen if the 
coupling interaction generated a pointer basis that var- 
ied randomly with time. This case corresponds to our 
3-noise simulation in which the noise fluctuates along 
all 3 directions simultaneously. Refs. and 01 show 
that a quantum phase transition occurs during QuAdS 
near t/T ~ 0.7. Before (after) this phase transition the 
quantum dynamics is essentially driven by Hi (H p ). As 
discussed in Section ITTTI1 the eigenstates of Hi (H p ) cor- 
respond to qubit spins aligned along x (z). 

• For noise with fluctuations along x and for t/T < 
0.7, [H(t), H int (t)] w and the instantaneous en- 
ergy eigenstates are essentially the same as the 
pointer basis states throughout the first 70% of 
the quantum evolution. Noise-induced dephasing is 
thus only effective during the final 30% of the evo- 
lution where [H(t), Hi nt (t)] ^ and the pointer 
basis differs from the instantaneous energy eigen- 
states. Only during the final 30% of the quantum 
evolution will noise act to dephase entanglement. 

• For noise fluctuations along y, [H(t) 1 Hi nt (t)] 7^ 
for all t, and the instantaneous energy eigenstates 
never correspond to the pointer basis. Thus noise- 
induced dephasing occurs throughout the entire dy- 
namical evolution. 

• For z-type noise, [H(t), H; mt {t)] w during the 
final 30% of the dynamical evolution, while noise- 
induced dephasing occurs over the initial 70% of the 
evolution. Thus dephasing of entanglement occurs 
over nearly 3/4's of the quantum evolution. 

These remarks suggest that z-type noise should have a 
larger impact on QuAdS than x-type noise, and that 
y-type noise should have the most severe impact of 
the three. As 3-type noise has fluctuations along all 3 
directions, it is unlikely that noise fluctuations of this 
type will remain aligned with any one of the coordinate 
axes throughout the entire adiabatic evolution. Thus 
3-type noise will sample the most damaging y-type 
fluctuations less often than y-type noise, and the more 
benign x-type fluctuations less often than x-type noise. 
As a result, one would expect 3-type noise to have a 
less severe impact on QuAdS than y-type noise, and 
more of an impact than x-type noise. This is what 
we see in our simulations (see Tables llT lVIf) . We see 
that knowledge of the occurrence of a quantum phase 
transition during QuAdS near t/T ~ 0.7, together with 
knowledge of how phase decoherence degrades pointer 
basis superpositions, allows us to understand most of 
the trends in our numerical results. Applying these 
arguments to a comparison of 3-type noise with z-type 
noise does not lead to any simple conclusion as far as 
we can tell. Our simulation results suggest that this 
case is in fact less straightforward as 3-type noise is less 
damaging to QuAdS performance than z-type noise at 
the lowest power levels simulated, but becomes more 
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detrimental to QuAdS at the larger power levels. 

(3) Using eq. (|54Jl , we can estimate the degree of decoher- 
ence in our simulations. From this formula we see that 
the least amount of decoherence will occur for a — a' = 1 
(i- e - d <t,<t' is largest): 

£>* ax = exp[-27 2 -rT^] . (55) 

For our simulations, a x — 0.2, r = 1, and 7 = 1. For 
x-noise with P = 0.005 and N = 13, < T(13) >~ 22 
and [H mt {t),H(t)] ^ for 0.7T < t < T, or 30% of the 
time for QuAdS to be carried out. Thus only 0.3T should 
be used in eq. I|55|l since only for this duration did the 
instantaneous eigenstates of H (t) differ from the pointer 
basis. Inserting these values into eq. (|53|l gives D^ax ~ 
0.590. Thus for our simulation of x-noise, QuAdS still 
possesses a substantial amount of quantum coherence. 
Repeating the above analysis for y-noise we again find 
eq. with a x — > a 2 .. Setting (a — a' — 1), we find: 

^L, = cxp[-2 7 2 TTa2] . (56) 

For P = 0.005, < T(12) >~ 35, and [H int (t), H(t)} ^ 
throughout all of the QuAdS. Thus, inserting T ~ 35 
into eq. Ij56(l . along with a y — 0.2, and the above 
values for t and 7 gives D-},^ ~ 0.06. Thus y-noise at 
P = 0.005 produces non-trivial decoherence effects. 

(4) Finally, we show that for weak noise D max provides 
an upper bound for the noise-averaged success probabil- 
ity of QuAdS. 

As explained in Section Hi Bl for a USA instance of EC3, 
the groundstate \E g (T)) of the final noiseless QuAdS 
Hamiltonian H(T) (eq. Q) encodes the unique solution 
to this instance. The probability P suc for QuAdS to suc- 
ceed on this instance is: 

P suc =TrP g p(T) , (57) 

where P g — \E g (T))(E g (T)\ is the projection opera- 
tor onto the final groundstate \E g (T)) and p(T) = 
\ip(T))(ijj(T)\ is the final density matrix. We can assess 
the impact of noise on the performance of QuAdS by 
evaluating the noise-averaged success probability P suc - 

P^ c = TrP g p(T) . (58) 

Since P g projects onto the final groundstate of the noise- 
less QuAdS Hamiltonian, it does not depend on noise 
and so: 

P^Tc = Tr P g p(T) 

= £ K(T)a;,(T)) {A^a^TKCT)}, (59) 

cr ,cr' 

where we: (i) have carried out the trace using the pointer 
basis states (eqs. and (J33J)); and (ii) used eq. and 



the definition of a CT (T) given below eq. I|44|) . Noting that 
the set of {a rT (T)a^ T , (T)} are the matrix elements of the 
density matrix p g (T) = \E g (T))(E g (T)\ in the pointer 
basis representation, we have 

Psuc = ^ Dcr'.cr (PgjT))^^ (PgjT))^ ;CT 

< D max Tr (p g (T)) 2 

< D max . (60) 

Here we have used that p g (T) describes a pure state so 
that {p g (T)) 2 = p g (T). From eqs. lO and (JHSl we see 
that the noise-averaged success probability for QuAdS is 
exponentially sensitive to the noise parameters appearing 
in D max . It is clearly of interest to extend the above 
analysis beyond the limit of weak noise. We leave this 
extension to future work. 



C. Future Work 

A number of directions for future work suggest them- 
selves. As mentioned above, we would like to extend 
our simulations to larger noise power levels to examine: 
(i) how rapidly the scaling exponent b increases with noise 
power; (ii) the degree to which the quality of the power- 
law and exponential fits continues to decrease with noise 
power; and (iii) how the quality of the exponential-fit 
compares with that of the power-law fit as noise power 
is increased. The aim of this latter point being to see 
if we can observe whether one of the two types of scal- 
ing laws begins to provide a significantly better fit than 
the other at some higher noise power level. We would 
also like to examine quantitatively how QuAdS perfor- 
mance is affected by noise which varies from one qubit- 
site to another. One would expect that, for this type 
of noise, noise-induced decoherence would be more effec- 
tive at hampering QuAdS performance than the uniform 
noise we have considered in this paper. It would also 
be interesting to extend the weak noise analysis of the 
upper bound on the noise-averaged success probability 
for QuAdS to stronger noise. Finally, Farhi et. al. [22| 
have argued that varying the path from Hi to H p might 
improve QuAdS performance. Noise can be thought to 
implement a random path variation, and so it might be 
interesting to examine whether a noise parameter regime 
exists where the performance enhancement of Ref. j2^ 
might occur. As pointed out in the Introduction, we did 
find noise realizations which did reduce the runtime of 
QuAdS on a given USA instance, but it appears that the 
predominant effect of noise at the power levels that we 
considered is to slow down QuAdS. It might be worth- 
while to examine whether noise-improved QuAdS might 
occur at small noise power levels. 
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VI. SUMMARY 

In this paper we have presented the results of a large- 
scale simulation of QuAdS in the presence of noise. We 
determined the noise-averaged median runtime (T(N)) 
for QuAdS to succeed in solving USA instances of the 
NP-Complete problem N-Bit Exact Cover 3. Clear evi- 
dence was found of the algorithm's sensitivity to noise. 
We simulated noise with 4 different types of polariza- 
tion, and our results are the outcome of approximately 
300,000 integrations of the Schrodinger equation. The 
scaling relation for (T(N)) versus the number of bits N 
was fit with both a power-law scaling (T(N)} = aN b and 
an exponential scaling (T(N)) = a [exp(feiV) — 1]. Both 
scaling relations provided excellent fits to the simulation 
results, although the quality of the fits were found to 
decrease slightly with increasing noise power. The qual- 
ity of the exponential fit decreased at a slightly slower 
rate than did that of the power-law fit. The variation of 
the scaling parameters a and b with mean noise power 
and noise polarization was determined. Our simulation 



results are summarized in Tables ITt lVIl These tables or- 
der the noise types according to which type most slowed 
down QuAdS. We also showed how noise- induced dephas- 
ing can cause decoherence in the dynamics of QuAdS, es- 
timate the amount of decoherence present in our simula- 
tions of nosiy QuAdS, and derive an upper bound for the 
noise-averaged QuAdS success probability in the weak 
noise limit that is appropriate for our simulations. 
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